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ABSTRACT 

We carried out ^^C0(1 - 0) and ^^C0(2 - 1) observations of 21 different regions in the vicinity of M86, NGC4438, and 
along the 120 kpc-long, Ha-emitting filamentary trail that connects them, aiming to test whether molecular gas can 
survive to be transferred from a spiral to an elliptical galaxy in Virgo's lO^K intracluster medium (ICM). We targeted 
Ha-emitting regions that could be associated with the interface between cold molecular clouds and the hot ionized 
ICM. The data, obtained with the 30m telescope of the Institut de Radioastronomie Millimetrique, led to the detection 
of molecular gas close to M86. CO gas with a recession velocity that is similar to that of the stars, — 265kms~^, and 
with a corresponding H2 mass of 2x10'' Mq, was detected ~10 kpc southeast of the nucleus of M86, near the peak of 
its HI emission. We argue that it is possible for this molecular gas either to have formed in situ from HI, or to have 
been stripped from NGC4438 directly in molecular form. In situ formation is nonetheless negligible for the 7x10® Mq 
of gas detected at 12:26:15.9+12:58:49, at ~10 kpc northeast of M86, where no (strong) HI emission is present. This 
detection provides evidence for the survival of molecular gas in filaments for timescales of ~100 Myr. An amount 
equivalent to 5x10^ A/0 of H2 gas that could be lost to the ICM or to neighboring galaxies was also discovered in 
the tidal tail northwest of NGC4438. A scenario in which gas was alternatively brought to M86 from NGC4388 was 
also examined but it was considered unlikely because of the non-detection of CO below or at the H I stream velocities, 
2000-2700 kms-^ 
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1. Introduction 

Stripping of gas during galaxy collisions can replenish the 
gas reservoirs of elliptical galaxies and reinitiate star for- 
mation in them, in particular in overdense cluster environ- 
ments. Cooling flows can further add to the deposition of 
external gas in cluster ellipticals (see Fabian, 1994, and 
references therein). The growth of these ellipticals is reg- 
ulated by the interplay of gas cooling and heating along 
gas filaments that are associated with either mergers or 
cooling flows and that have often been discovered in the 
Ha and [Nil] images of local clusters (Cowie et al., 1983; 
Heckman et al., 1989; Conselice et al., 2001; Kenney et al., 
2008). This interplay manifests itself and can be examined 
through the detection of multi-phase gas emission lines in 
the same clouds. Indeed, warm atomic gas, warm molec- 
ular gas, and cold molecular gas clouds have been found 
to coexist in spatially resolved filaments around NGC1275 
(Perseus A), where they were detected through Ha, H2, 
and CO emission (Conselice et al., 2001; Salome et al., 2006, 
2011; Johnstone et al., 2007). 

Several processes, however, occur in cluster environ- 
ments that can impede the transfer of gas from one galaxy 
to another. Ram pressure stripping (Gunn & Gott, 1972), 
for example, can disperse the diffuse atomic gas into the 
intracluster medium (ICM). Even though the ram pressure 
drag is negligible for the motions of the dense molecular gas 
clouds (Nulsen, 1982; Kenney & Young, 1989), it can even- 



tually destroy them by depleting the atomic gas reservoir 
from which they reform (Growl et al., 2005). The molecu- 
lar clouds are also several orders of magnitude cooler than 
their environment and can be easily destroyed by X-rays 
if they are not self-shielded. The ICM X-ray spectrum of- 
ten peaks in the range 0.7-2.0 keV, indicating temperatures 
of ~10''K (e.g., Rangarajan et al., 1995; Machacek et al., 
2004; Fabian et al., 2006; Randall et al., 2008; Tamura et 
al., 2009). 

Evidence that molecular clouds do nonetheless survive 
or reform in >10^K intracluster and intragalactic media has 
been found. Braine et al. (2000) discovered GO in a tidal 
tail dwarf galaxy of Arpl05, in the X-ray emitting medium 
of Abell 1185 (Mahdavi et al., 1996). Because the CO- 
emitting regions spatially coincided and had comparable 
kinematic properties with the Hl-emitting regions, Braine 
et al. (2000) proposed in situ formation of molecules. In 
addition to the molecular gas detected in the cooling-flow- 
related filaments around NGC1275 in Perseus (Salome et 
al., 2006, 2011; Johnstone et al., 2007), Vollmer et al. (2008) 
found extraplanar CO near NGG4522 in Virgo, in regions 
where atomic gas was brought by ram pressure. In ES0137- 
001, located in Abell 3627, a ~10'^K ram-pressure-driven 
gas tail has active star formation and H2 emission (Sun et 
al., 2010; Sivanandam et al., 2010). In the compact group 
of Stephan's Quintet, H2 molecules not only survive in a 
medium of comparable temperature (Trinchieri et al., 2005; 
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Fig. 1. Finding chart of our observations. The positions of previous CO observations (Wiklind et al. 1995; Braine et al. 1997) 
are marked with light green circles. HI line contours (dark green) and X-ray contours (brown) are plotted over the narrow band 
Ha-|- [Nil] image of Virgo's core (adapted from Kenney et al. 2008). The velocities appearing in gray correspond to those of the 
Hq gas as measured by the same authors. In red, we mark the Ha clouds that we observed searching for associated CO emission. 
Upon detection, we also indicate the velocity of the molecular gas traced by the ^^^^00(1 — 0) line. In addition to these 20 regions, 
we also observed at 12:27:30.0+13:12:12.0, northwest of NGC4435, in the clouds that Cortese et al. (2010a) denote as plume. 



O ' Sullivan et al., 2009), but they are mainly responsible for 
the gas cooling instead of the X-rays (Cluver et al., 2010). 

An excellent system to test whether molecular gas could 
be directly brought to an elliptical galaxy, under extreme 
conditions, from spiral(s) interacting with it is in the Virgo 
cluster. The dust clumps that are detected in the giant el- 
liptical galaxy M86 (also known as NGC4406) could have 
originated from its past interaction with its neighboring 
spiral galaxy NGC4438 (Cortese et al., 2010b; Gomez et 
al., 2010). Dust and gas could have been captured through 
ram pressure stripping of the interstellar medium (ISM) of 
NGC4438 by that of M86 (Kenney et al., 2008), through 
dynamical instabilities, or through both processes. The HI 
content deficit of NGC4438 relative to its own CO emission 
confirms that ram pressure has acted on its ISM (Vollmer et 
al., 2005, 2009; Hota et al, 2007). Dynamical perturbations 
continue to distort NGC4438, as indicated by the presence 
of stars and gas in its northern tidal tail toward NGC4435 
(Combes et al., 1988). The interaction between NGC4438 
and M86 is confirmed by a spectacular, ~120kpc-long, Ha- 
emitting tidal bridge that connects them (Kenney et al., 
2008). The smooth Ha velocity gradient along this bridge 
convincingly places it to the ICM of Virgo, which is rich 
in X-ray emission (Forman et al., 1979; Rangarajan et al., 
1995; Machacek et al., 2004; Liu & Bregman, 2005; RandaU 
et al., 2008). Since the X-ray heating and the ram pressure 
disperse the ICM atomic gas, a significant fraction of the 
Ha emission could originate from the ionized gas in the 



outer layers of molecular clouds (Ferland et al., 2009). This 
makes the Ha clumps in the NGC4438-M86 bridge prime 
targets for testing whether star formation can be initiated 
from gas that has either been brought from NGC4438 to 
M86 in the molecular state, or from molecules that recom- 
bined in situ in M86 after being transferred in the atomic 
state. 

We obtained mm data with the 30 m telescope of the 
Institut de Radioastronomie Millimetriquc (IRAM) to test 
for the presence of molecular gas in some of the Ha fila- 
ments that are embedded in Virgo's hot ICM. We focused 
on Ha emission peaks, aiming to determine whether the 
warm atomic gas emission could be associated with cold 
gas clouds, to compute the amount of molecular gas and 
potential star formation rate (SFR) in M86, and to test the 
gas origin by examining whether the molecules could have 
been transferred there or formed in situ. We assume a dis- 
tance of 17.5 Mpc to Virgo throughout this work (Mei et 
al., 2007). 

2. Observations and data analysis 

The observations were carried out with the 30m telescope 
at Pico Veleta, Spain. We observed twenty Ha-emitting re- 
gions in or near NGC4438 and M86, and along the tidal 
bridge that connects them (Figure 1), to test whether they 
contain detectable amounts of CO. In addition to these po- 
sitions, we also observed clouds in the diffuse. X-ray emit- 
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Fig. 2. Spectra of tlie plume and of the Virgo regions shown in Figure 1 with ^^CO(l — 0) line detections. Gaussian fits to the CO 
line profiles (black) are overplotted in red. 
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Table 1. ^^CO(l — 0) and ^^C0(2 — 1) line properties and derived gas masses in different Virgo regions. 



Position 


Coordinates 


Integration 


Line 


5(CO)» 




FWHM° 


^^ d 

-tmb 


resolution 






(.72000) 


time (min) 




(.Jy km s^^) 


(kms^^) 


(kms^^) 


(mK) 


(kms^^) 


(Me) 


NGC4438 


12:27:15.60+13:00:32.0 


15 


"^CO(l-O) 


226±6 


72+6 


395+12 


108+14 


10 


5.7x10'* 


NGC4438 




15 


1^00(2-1) 


208±8 


61+7 


362+14 


108+15 


13 




Al 


12:27:43.70+12:59:22.8 


110 


^^CO(l-O) 


3.5+0.6 


-76+5 


46+11 


13.4+3 


21 


8.8x10" 


Al 




109 


^^C0(2-l) 








<5.1 


47 




A2 


12:27:42.10+12:58:51.4 


25 


^^CO(l-O) 


5.4+0.6 


-4+2 


37+5 


27+5 


10 


1.4x10'' 


A2 




25 


^^C0(2-l) 


3.8+0.6 


-12+3 


29+6 


24+5 


10 




El 


12:27:38.60+12:59:52.0 


90 


^^CO(l-O) 


2.3+0.4 


6+5 


52+12 


7.9+1.6 


21 


5.8x10*' 


El 




90 


^^C0(2-l) 








<3.8 


47 




EIW 


12:27:37.40+12:59:38.3 


37 


^^CO(l-O) 








<5.4 


47 


<3.2xl0'' 


EIW 




37 


i^CO(2-l) 








<3.9 


47 




E2 


12:27:44.00+13:03:00.0 


50 


^^CO(l-O) 


8.8+0.7 


144+2 


42+5 


39+3 


10 


2.2x10'' 


E2 




50 


i^CO(2-l) 


3.8+0.7 


141+2 


21+6 


29+5 


13 




E2E 


12:27:45.50+13:02:48.0 


65 


^^CO(l-O) 


8.7+0.6 


176+1 


38+3 


42 ±4 


10 


2.2x10'' 


E2E 




65 


i^CO(2-l) 


5.3+0.6 


172+2 


28+5 


33+4 


10 




E2W 


12:27:42.10+13:03:15.0 


80 


^^CO(l-O) 


2.6+0.6 


157+10 


84+22 


6.1+2.1 


21 


6.5x10" 


E2W 




80 


^^C0(2-l) 








<4.8 


47 




V5 


12:27:31.60+13:01:38.0 


42 


^^CO(l-O) 








<5.4 


47 


<3.2xl0" 


V5 




42 


i^C0(2-l) 








<3.8 


47 




VI 


12:26:56.20+12:59:40.0 


52 


^^CO(l-O) 








<4.7 


47 


<2.8xl0" 


VI 




52 


^^C0(2-l) 








<6.4 


47 




V2 


12:26:19.00+12:58:10.0 


43 


^^CO(l-O) 








<4.7 


47 


<2.8xl0" 


V2 




43 


1^00(2-1) 








<4.1 


47 




V7 


12:26:15.90+12:58:49.0 


50 


^^CO(l-O) 


2.8+0.5 


-527+6 


59+12 


8.4+1.9 


21 


7.1x10" 


V7 




45 


i^CO(2-l) 


2.6+0.4 


-495+6 


66+13 


7.0+1.4 


21 




V4 


12:25:47.50+13:02:44.0 


23 


^^CO(l-O) 








<6.5 


47 


<3.9xl0" 


V4 




23 


i^CO(2-l) 








<4.4 


47 




V3 


12:26:17.00+12:54:54.0 


37 


^^CO(l-O) 








<5.3 


47 


<3.1xl0" 


V3 




37 


i^CO(2-l) 








<5.6 


47 




V3W 


12:26:14.90+12:54:54.0 


129 


^^CO(l-O) 


4.9+0.6 


-265+6 


89+15 


10+2 


21 


1.2x10'' 


V3W 




129 


i^CO(2-l) 








<4.4 


47 




V3WW 


12:26:12.80+12:54:54.0 


115 


^^CO(l-O) 


3.2+0.4 


-265+4 


57+9 


11+2.3 


10 


8.1x10" 


V3WW 




115 


^^CO(2-l) 


3.1+0.5 


-263+6 


68+16 


8.5+2.2 


13 




V3NW 


12:26:15.50+12:55:22.0 


40 


^^CO(l-O) 








<8.0 


47 


<4.7xl0" 


V3NW 




40 


1^00(2-1) 








<5.5 


47 




V6 


12:26:15.50+12:55.50.0 


35 


^^CO(l-O) 








<7.5 


47 


<4.4xl0" 


V6 




35 


^^C0(2-l) 








<4.5 


47 




V8 


12:26:08.00+12:56:12.0 


35 


^^CO(l-O) 








<7.3 


47 


<4.3xl0" 


V8 




35 


1^00(2-1) 








<6.5 


47 




M86SW 


12:26:11.40+12:56:15.0 


23 


^^CO(l-O) 








<12 


47 


<7.1xl0" 


M86SW 




23 


1^00(2-1) 








<12 


47 




PLUME 


12:27:.30. 00+13:12:12.0 


10 


i2(X)(l-0) 


8.1+0.2 


-14 


1.2 


1018+27 


0.8 




PLi::iIE 




10 


'-('()(2-l) 




-il 


i.2 


782=27 


0.1 





Fluxes obtained from the Gaussian fitting to the hne profiles, and multiplied by a temporature-to-flux conversion factor of 
5.0 Jy/K that is appropriate for the 30m telescope at all observed frequencies. '''' All velocities are heliocentric. ^'^^ Resolution- 
corrected line width. All limits are computed at 3a levels. ^""^ A standard CO-to-H2 Galactic value of 4.6 Mq (K kms^^pc^)"^ 
was used for the mass computation. The 3a limits were computed for a width of 47 km s~ ^ . This is close to the median CO (1 — 0) 
width for all regions but the plume and NGC4438. 



ting plume northwest of NGC4438 and NGC4435. which 
is most likely attributed to Galactic cirrus (Cortese et al., 
2010a). 

We simultaneously observed the -"^^CO J=l— and 
3—2—1 lines with the EMIR receivers. The receivers were 
tuned to the frequencies of the two transitions at 115.271 
GHz and 230.538 GHz, respectively, shifted to the velocity 
of the Ha gas in each position. All of the 4MHz, 'WILMA, 
and VESPA backcnds, with corresponding resolutions of 4 
MHz, 2 MHz, and 0.3 MHz, were simultaneously attached 
to allow consistency checks, with the WILMA data being 
primarily used in this work. VESPA data are only shown 
for the plume. The observations were performed in wobbler- 
switching mode with a throw of 2.'40. The telescope point- 
ing was monitored by observing a bright radio source every 
two hours, indicating a pointing accuracy of 2" to 3". The 



system temperature varied from 270 to 500 K at 3 mm. We 
typically spent 0.25-2 hours on source, which led to la noise 
levels of 2 mK to 10 mK at the frequency of the ^^C0(1 - 0) 
transition (Table 1). 

The data were reduced with the CLASS software of 
IRAM. Wc! avciragcid the spectra of each region after re- 
moving bad channels and subtracting linear baselines, drop- 
ping scans with baseline fluctuations. The resulting spectra, 
binned to 10-21 kms""'^ channels, are shown in Figure 2. 
For the plume, the high-resolution VESPA spectra are un- 
binned. A Gaussian function was fitted to each line for 
the derivation of its width, peak position, temperature T, 
and flux S, which are all shown in Table 1. Assuming that 
the observed clouds are unresolved, T corresponds to the 
main beam temperature. It was converted to a flux using a 
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Fig. 3. Our CO observation positions (red circles) and positions of previous CO observations (Wiklind et al. 1995; Braine et al. 
1997; green circles), over Herschel 250 ^m SPIRE imaging data (adapted from Cortese et al. 2010b). The Ha intensity contours 
are overplotted in blue. 



temperature-to-flux conversion factor of 5.0 Jy/K (appro- 
priate for both the 1 and 3 mm observations). 

The cold H2 gas masses, which are also presented in 
Table 1, were computed as 



Ml 



H2 



23.5 /co ^bD 

(1 + Z)3 



M, 



(1) 



where /co is the 115GHz line intensity in K kms^^, fis 
is the main beam area in arcsec^ (for unresolved clouds), 
-Dl is the luminosity distance in Mpc, and a is the CO 
luminosity to H2 mass conversion factor (Solomon et al., 
1997). We used a standard Galactic conversion factor of 
4.6M0 (K kms-ipc^)-^ (Bloemen et al., 1986; Downes 
et al., 1993), and a telescope beam of 22" (i.e., 1.8 kpc) 
at 3 mm. The 3cr mass limits for non-detections were 
computed using a resolution (and width) of 47 kms~^, 
which is characteristic of the clouds under examination (see 
Table 1). 



3. Results: CO-emitting vs non-emitting regions 

Narrow ^"^CO^i - 0) and 1^00(2 - 1) emission is found in 
the plume clouds northwest of NGC4438 and NGC4435, at 
12:27:30+13:12:12. Both line widths are of only 1.2 kms"\ 
confirming that the emission originates from Galactic cir- 
rus clouds. This possibility was suggested by Cortese et al. 
(2010a) based on the ^^C0(1 — 0) line width at a nearby 
position, 12:27:30+13:12:29. CO lines were detected in an- 
other 10 of the 20 regions shown in Figure 1, but their 
kinematics were atypical of cirrus clouds, with line width(s) 
exceeding 30 kms~^ in all cases. A cross examination of the 
high-resolution VESPA data confirms this finding. 

Several regions in Virgo, at ~10 kpc from the _B-band 
center of M86 (at 12:26:12+12:56:49) had a CO detection. 



The region V7, which is located 11 kpc northeast of M86 
and in the northwest tip of the local Ha filamentary struc- 
ture (Fig. 1), contains 7.1xlO^M0 of gas that is moving at 
-527kms~^ This is very close to the measured Ha reces- 
sion velocity, — 500kms~^ (Kenney et al., 2008). This line- 
of-sight velocity precludes a clear answer to the question of 
whether the clouds are bound to M86. The escape velocity 
at a distance of 10 kpc from the center of M86 is high, 
470 kms^^, because its bulge mass equals 2.6x10^^ Mq. 
The latter is given in units of 10^° Mq by C2crioo^-Reff , where 
C2 is a unitless coefficient depending on the matter distri- 
bution, (Tioo is the stellar velocity dispersion in units of 
100 km s^^, and i?off is the effective radius in kpc. The value 
of C2 for giant ellipticals is 1 (Bender et al., 1992). For M86, 
(Tioo is 2.2 (Smith et al., 2000) and i?cff is 5.5 kpc (Gavazzi 
et al., 2005) in the B band. 

CO gas was also found in the V3W and V3WW re- 
gions, along an Ha stream that results in the center of 
M86, at 9 to 10 kpc away from it (Figure 1). The CO emis- 
sion coincides with the H I emission peak of M86 (Li & van 
Gorkom, 2001; Kenney et al., 2008). The recession veloc- 
ity of the CO, — 265kms~"'^, agrees well with that measured 
from the HI, — 240kms~^ (Bregman & Roberts, 1990), and 
from the Mg stellar absorption features in optical wave- 
lengths, -244kms-i (Smith et al., 2000). The i2C0(l - 0) 
FWHM in these regions, 50-90 km s~^, is comparable with 
the narrow HI component FHWM, eOkms"^ (Bre gman & 
Roberts, 1990). The combined H2 mass for both V3W and 
V3WW regions is 2.OxlO''M0. If stars are forming out of 
this gas, and if the galaxy-integrated SFR calibration of 
Gao & Solomon (2004) is applicable for these clouds, then 
the SFR will be 0.03 M0yr^^. Likewise, the combination 
of the modified black body fit to SPIRE data (Cortese et 
al., 2010b, Figs. 3 and 4) of the same regions that was per- 
formed by Gomez et al. (2010) and the Kennicutt (1998) 
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relation lead to a dust-based SFR estimate that is also on 
the order of O.OlMoyr"^. 

No molecular gas was found at their neighboring V3, 
V3NW, V6, V8, and M86SW regions. The latter is along 
the same Ha stream connecting V3W and V3WW to the 
very center of M86. For T<12mK and for a line width 
of 47kms-\ the gas mass in M86SW is <7x106Mq. 
Molecular gas has not been detected at the center of M86 
either (Wiklind et al., 1995; Braine et al., 1997); if it exists, 
its mass must be below 5x10^ Mq for the assumptions in 
this work. 

No gas was found either at VI, V2, and V5 of the 
tidal bridge between NGC4438 and M86, or at V4, which 
is located in the figure-of-eight shaped Ha loop north- 
west of M86^. If gas does exist in any of these regions, 
its mass is below 4x10^Mq. Previous observations by 
Wiklind et al. (1995) at another location of this loop, 
12:26:02-1-12:59:21, also led to a non-detection, with an up- 
per limit of 7.9x10'^Mq. 

Most of the CO gas is detected in or around NGC4438. 
At the center of NGC4438, the telescope beam encom- 
passes 6xlO*M0 of H2 gas that has a recession velocity 
identical to that of the Ha. The 1^00(1 - 0) line profile 
also agrees well with that found from previous observations 
by VoUmer et al. (2005). Both the 1^00(1-0) and the 
""^^00(2 — 1) profiles indicate that the gas is in a disk whose 
rotational velocity reaches 200 km s^^ within the inner 4 
kpc of NGC4438^. New observations are presented here for 
the gas clouds at El, and in the regions E2, E2E, and E2W 
at ^10 kpc north of the nucleus. We specifically observed 
parts of the Ha filament that is oriented perpendicularly to 
the tidal tail northeast of NGC4438, which was previously 
observed by VoUmer et al. (2005). The nearest positions 
of our observations and those of VoUmer et al. (2005) are 
one beam (22") away. The clouds in the E2 complexes are 
moving at -1-70 to -1-100 km s""'^ from systemic velocity along 
the line of sight. Their total H2 mass is 5.1x10'^Mq. Their 
counterparts in the south of the nucleus are named Al and 
A2. They contain 2.3x10^ Mq of H2 gas, and they are mov- 
ing at —80 to —150 kms^^ from systemic velocity. For both 
the northern and southern clouds, the observed velocities 
are consistent with gas that is either part of a rotating disk, 
or experiencing tidal and ram pressure stripping (VoUmer 
et al., 2005). If all of these A and E2 cloud complexes are 
stripped, a considerable (14%) fraction of the gas in the 
center of NGC4438 (see Table 1) is potentially lost to the 
ICM and to neighboring galaxies. 

4. Discussion: the origin of the molecular gas near 
M86 

Li & van Gorkom (2001) attributed the origin of the HI 
emission, which peaks near CO-emitting regions, to gas 
that cooled down from a compressed, shocked, and heated 

^ This loop belongs to Virgo's ICM and part of it could be 
bound to M86, because its Ha recession velocity ranges from 
-230 to -450 kms"V 

^ The differences in the observed ^^C0(1 - 0) and^2CO(2 - 1) 
line properties of other regions, such as E2W, E2, and E2E, can 
be accounted for by differences in the gas density and excitation 
properties, or differences in the sensitivity and the beam size of 
the observations at the two frequencies, if some filaments are 
resolved. 




Fig. 4. Our observations plotted over a 4'x4' zoom around M86 
at 350 ^m. The far-infrared image is taken from the first public 
data release of the Herschel Virgo Cluster Survey (HeVICS), 
which was initially presented in Davies et al. (2010). The CO 
positions are color-coded as in Figure 3. 



medium. In cooling flows, however, the ICM gas is depleted 
of its dust due to sputtering by X-rays. Dust sputtering 
leads to gas-to-dust ratios that exceed their typical values 
for local galaxies (i.e. several hundreds; Young & Scoville, 
1991; Wikhnd et al, 1995) in timescales of only Myr in 10*^- 
lO^K environments (Clemens et al., 2010). On the other 
hand, the presence of dust in all M86 regions where CO 
was detected (Figs. 3 and 4; Cortese et al., 2010b; Gomez 
et al., 2010) indicates that the gas and the dust were com- 
monly accumulated. Their relative mass ratio indeed favors 
the capture from an external galaxy rather than the forma- 
tion from a cooling ICM. For the inner 2' of M86, the total 
15-20 K dust mass is {2-5)xl0'^ Mq (Gomez et al., 2010). 
If the molecular gas distribution follows the dust distribu- 
tion, most of the H2 gas in the inner 2' will be accounted 
for by the observations at V6, V3NW, V3W, V3WW, V8, 
M86SW, and at the center of M86. The total H2 mass 
will then be in the range (2-5)xlO^M0. Furthermore as- 
suming that the center of M86 comprises 1.1x10^ Mq of 
atomic gas based on its Ha-|- [N II] luminosity (Kenney et 
al., 2008), and another 5.6 xlO''M0 based on its HI ob- 
servations (Bregman & Roberts, 1990), we deduct that its 
gas-to-dust ratio will be in the range 20-60. The range can 
be broader when considering the different beam sizes of the 
observations. For comparison, Cortese et al. (2010b) give a 
10-20 K dust mass of 2x10^-2x10'^Mq, leading to a gas- 
to-dust ratio of 30-300 (including helium) for a dust knot 5 
kpc northwest of NGC4438. 

With its Ha-emitting traU, NGC4438 is the best can- 
didate for the external capture of the gas and the dust in 
M86. While the trail started elongating as the two galax- 
ies were receding from each other, most of the gas and the 
dust was likely captured when the galaxies where near the 
pericenter of their encounter. If gas was transferred from 
NGC4438 to M86 in ionized or neutral state through dy- 
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namical instabilities or ISM-ISM stripping, CO could be 
detected in M86 if it reformed in situ. Molecular forma- 
tion is possible to occur in situ in the V3W and V3WW 
complexes. The column density threshold for the efficient 
conversion of HI to H2 is in the range 3xl0^°-10^^ cm~^ 
(Schaye, 2004). The Hi-deduced column density in these 
regions is close to the 1.6 x lO^^cm"^ contour of Kenney 
et al. (2008). Likewise, the CO-based H2 column density 
is 2.9x10^0 and LOxlO^Ocm^^ for V3W and V3WW, re- 
spectively, for the adopted CO-to-H2 Galactic conversion 
factor. The HI column density requirement would be satis- 
fied if the HI-to-H2 mass ratio was >3, which is comparable 
to that of both moderately stripped spirals and ellipticals 
(Kenney & Young, 1989; Wiklind et al., 1995). The present 
computations assume that the gas emission is smeared out 
across each telescope's beam. When taking into account 
the clumpy distribution of the clouds, the column density 
requirement can be easily met. For this scenario to hold, 
the CO needs to be detected in regions where dust is also 
present, because the dust particles act like catalysts for the 
HI atoms to recombine upon. This condition is also satis- 
fied. 

Contrarily to the V3W(W) complexes, the bulk of the 
molecules could not have formed in situ in V7, where the 
Hl-deduced column density is outside the last contour of 



Kenney et al. (2008), i.e., below lO^^cm^^. The detection 
of CO there provides evidence for the survival and trans- 
port of self-shielded molecular clouds in the hot ICM of 
Virgo. The molecules could have started forming, e.g., near 
the H I peak south of M86 and then rotated around it in a 
half-circular orbit of ^30 kpc. If the clouds were moving at 
300 km s^^ (i.e. close to the V3W and V7 difference in reces- 
sion velocity), their travel time would be 100 Myr. Because 
the formation of new stars from H2 in isolated clouds takes 
place in timescales of 10^-10^ yr (e.g., Krumholtz et al., 
2011), this result implies that the collapse, the fragmenta- 
tion, and the reformation of the molecular clouds from gas 
in their vicinity is a continuous process inside filaments. 

Alternatively, gas clumps could have been transferred 
from NGC4438 to MSG directly in a dense, molecular state. 
The molecular gas has then survived near M86 for a time 
comparable to the travel time of NGC4438 from the peri- 
center to its present location. The orbital simulations of 
VoUmer et al. (2009) also suggest that NGC4438 has been 
moving through the hot X-ray gas for over only 100 Myr 
after its encounter with MSG. This timescale is insufficient 
for the various energy /momentum transport processes that 
act on the molecular clouds to fully destroy them. These 
processes include the heat conduction from the surround- 
ing medium, the viscous flow stripping, and the Kelvin- 



7 



Dasyra et al.: CO in Virgo's hot intracluster medium 



stack of V6+V3NW+V3+V3W+V3WW+V8+M86SW 



CO(l-O) 

0^004 h CO(l-0) + (2-l) 



0.002 - 



0.002 




-2000 -1000 1000 

Velocity (kms~^) 



2000 



3000 



Fig. 6. Stacked spectra of all positions south of M86. No gas 
is detected at any velocity other than the galaxy's recession 
velocity. 



Helmholtz instabilities, with the heat conduction leading 
to the highest mass loss rates in this environment (Nulsen, 
1982). In the ICM, the mean free path of the electrons 
is large compared to the radii of the molecular clouds, 
leading to saturated heat conduction and long evaporation 
timescales (Cowie & McKee, 1977). For a cloud with ra- 
dius of 10 pc and density n of 10'^ cm~'^ that is embedded 
in a 10 ''K medium with density of 10""* cm~'^, the evapora- 
tion time is 3 Gyr. This time can increase when magnetic 
fields and radiative cooling are taken into account (Cowie 
& McKee, 1977; McKee & Cowie, 1977), because cooling 
is very efficient above lO^K (Dalgarno & McCray, 1972). 
It effectively dominates over heating for clouds that are 
stripped outside of galaxies and that are not sufficiently 
heated by stellar ultraviolet radiation (VoUmer et al., 2001). 
The atomic gas will then cool down, condense, and reform 
II2 within ~10^/n yrcm"^ (HoUenbach & Tielens, 1997). 
It is therefore plausible that the molecular gas observed in 
M86 was directly stripped from NGC4438, even though it 
was not detected in our observations of, e.g., the VI region. 
This might occur if there is no strong source of gas excita- 
tion, such as star formation, shocks, cosmic rays (Ferland et 
al., 2008), or diffuse intracluster light (Mihos et al., 2005). 

An entirely different scenario that remains to be tested 
is whether the gas in M86 could originate from NGC4388, a 
spiral galaxy at ~100kpc in projection on its south, with an 
estimated II2 reservoir of 10^-10^ Mq (Kenney & Young, 
1989; Papadopoulos k Seaquist, 1998). NGC4388 could be 
interacting with the M86-NGC4438 complex, if all galax- 
ies were at a sufficiently close distance (e.g., Willick et al., 
1997). The interaction scenario was considered after the 
detection of a long (>120 kpc) trail of neutral gas starting 
from NGC4388 and heading in the direction of M86 and 
NGC4438 (Oosterloo & van Gorkom, 2005, Fig. 5), and 
by the strong gas stripping in this system (e.g., Kenney 
& Young, 1989; VoUmer et al., 2009). Very Large Array 
data indicated that NGC4388 has a radio jet that is driven 
by an active nucleus (Colbert et al., 1996; Irwin et al., 
2000) and that could be adding to the gas expulsion mech- 
anisms. A massive outflow has been observed in its Ha 
and [Nil] images (Veilleux et al., 1999; Kenney et al., 
2008), with an Ha velocity ranging between 2000 kms""'^ 
and 2400kms~i (Yoshida et al., 2004). The velocity of the 
neutral, H I-emitting gas along the trail is likewise between 



2000-2700 km s-i (Cayatte et al, 1990; Oosterloo & van 
Gorkom, 2005; Chung et al., 2009). The projected H I detec- 
tion limits are close to the Ha detection limits in the vicin- 
ity of M86 (Fig. 5). However, the non-detection of CO mov- 
ing at any velocity between — 250kms~^ and 2700 km s~^ 
in the south of M86 (either in individual or stacked spectra; 
Fig. 6) renders this scenario improbable. A similar conclu- 
sion was drawn by Oosterloo & van Gorkom (2005) based 
on the HI gas kinematics. Additional (and deeper) observa- 
tions would be needed to test the viability of this scenario. 

5. Summary 

Using IRAM's 30m telescope, we investigated the presence 
of molecular gas in Virgo's hot ICM, and specifically in the 
vicinity of M86, NGC4438, and in the 120 kpc-long tidal 
trail that connects them. We focused on Ha emission peaks, 
which could originate from the interface of dense clouds 
with self-shielded CO molecules and the X-ray emitting 
ICM, because ram pressure stripping has dispersed most 
of the diffuse ionized gas. We observed 21 regions in the 
i2C0(l - 0) and 1^00(2 - 1) transitions and we draw the 
following main conclusions. 

— Narrow CO emission, with line widths of only 1.2 
kms"\ was detected at 12:27:30+13:12:12, northwest 
of NGC4438 and NGC4435. It is tracing Galactic cirrus 
clouds. In no other position were the CO line widths 
characteristic of those in the cirrus. 

— CO was detected at more than 10 kpc away from the 
centers of NGC4438 and M86. For example, strong 
CO emission was found in the tidal structure north of 
NGC4438, in the regions E2E, E2, and E2W. Through 
these cloud complexes, 5x10''Mq or 9% of the gas in 
the center of NGC4438 could be lost to the ICM and 
potentially to neighboring galaxies. 

— CO was discovered in an Ha stream that ends in the 
center of M86, at ^^10 kpc southeast from the center, 
and near the peak of the H I emission. The CO recession 
velocity agrees well with that of the H I and the stars. 
The ^^C0(1 — 0) line properties indicate the presence 
of 2x10^ Mq of H2 in the combined regions V3W and 
V3WW. Some of this gas could have formed in situ given 
that the local HI surface brightness exceeds 10^" cm~^. 

— The filamentary structure 11 kpc northeast from the 
center of M86, in the region that is denoted V7, con- 
tains CO moving at — 530kms^^. Its corresponding H2 
mass is 7x10^ Mq. Because there is no associated HI 
emission in this region or in its vicinity, we attribute 
this CO detection to gas that was primarily transferred 
to V7 in molecular form. We interpret it as survival of 
self-shielded, re-forming molecular clouds in filaments 
for ^100 Myr. This is longer than the fragmentation 
timescale of individual isolated clouds, and shorter than 
the evaporation timescale due to heat conduction from 
the ICM. 

— The integrated gas-to-dust ratio in the inner 2' of M86, 
20-60, is too low to be consistent with gas cooling from 
the ICM. The molecular gas detected in and near M86 
is therefore likely to have originated from an external 
galaxy. 

— The most probable scenario for the molecular gas origin 
of M86 is that it comes from NGC4438. The best evi- 
dence for this remains the smooth Ha velocity gradient 
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in the bridge between them, combined with the similar 
Ha and CO velocities, where detected. 
— The scenario of gas being transferred to M86 from 
NGC4388 on its south was also examined, given the H I 
stream that is seen in projection between these galax- 
ies. However, the non-detection of CO below or in the 
observed range of HI recession velocities, i.e., between 
2000—2700 kms~^, does not favor this scenario. 
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